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We study the dynamics of the entanglement between two oscillators that are initially prepared 
in a general two-mode Gaussian state and evolve while coupled to the same environment. In a 
previous paper we showed that there are three qualitatively different dynamical phases for the 
entanglement in the long time limit: sudden death, sudden death and revival and no-sudden death 
[Paz & Roncaglia, Phys. Rev. Lett. 100, 220401 (2008)]. Here we generalize and extend those 
results along several directions: We analyze the fate of entanglement for an environment with 
a general spectral density providing a complete characterization of the evolution for ohmic, sub- 
ohmic and super-ohmic environments. We consider two different models for the interaction between 
the system and the environment (one where the coupling is through position and another where 
the coupling is symmetric in position and momentum). Finally, we show that for non-resonant 
oscillators the final entanglement is independent of the initial state and that it may be non-zero at 
very low temperatures. 



I. INTRODUCTION 



The creation and manipulation of entanglement is an 
important issue not only because of its fundamental im- 
plications but also due to its practical applications. In 
fact, during the last decade many new discoveries re- 
garding the physics of entanglement were made [T] . The 
development of quantum algorithms and cryptographic 
schemes were probably the driving forces behind most 
of the research on entanglement. In fact, nowadays en- 
tanglement is regarded not only as a peculiar feature of 
quantum systems but also as a physical resource. Entan- 
glement manipulation was studied first for finite dimen- 
sional systems but later continuous variable systems [2] 
were considered. In fact, several experiments showed the 
successful implementation of quantum teleportation [5] 
and cryptographic protocols [4] for such systems. In this 
context, it is important to take into account the effects in- 
duced by the interaction between a composite (eventually 
entangled) quantum system and its surrounding environ- 
ment. Indeed, decoherence in some cases can be devas- 
tating: Thus, due to the interaction with the environ- 
ment, entanglement within a composite system can dis- 
appear in a finite time. This phenomenon, that was first 
discussed and analyzed for systems made out of qubits 
[SJ ini E] became known as "sudden death" of entangle- 
ment (SD). But the fate of entanglement for a quantum 
open system is not at all evident and some surprising re- 
sults were also obtained: For example, it was shown that 
under certain conditions the environment can act as a 
quantum channel thorough which entanglement can be 
created |S]. In this case, even if the initial state of the 
system is separable, the final state could be entangled. 
A large number of recent papers study this and other is- 
sues that characterize the dynamics of entanglement in 
systems of qubits interacting with common or indepen- 
dent environments [H El UHl El US HH [14] . On the other 
hand, continuous variable systems were also investigated 



and similar results emerged. For example, the degrada- 
tion of entanglement for harmonic systems interacting 
with bosonic reservoirs was analyzed [151 [Ml [17] . Also, 
the fate of initially entangled states (two-mode squeezed 
states) interacting with a common bath was studied using 
different approximations [TS] [13 [20] ■ Among the inter- 
esting results that emerged from those works it is worth 
mentioning that in |19j a condition for the existence of 
sudden death was deduced under the RW-Markovian as- 
sumption. More recently, the non-Markovian regime was 
also analyzed [211 [221 [23] and a zoo with different long 
time behavior emerged. Thus, it was realized in [21 that 
non-Markovian effects modify the condition obtained in 
[19] for the existence of SD. 

More recently, [21], we provided a unified picture of 
the different qualitative dynamics of entanglement for 
general gaussian states in non-Markovian environments. 
There, we showed that the asymptotic dynamics of en- 
tanglement can be described by three possible phases: 
SD (sudden death) , SDR (sudden death and revival) and 
NSD (no sudden death). The existence of an exact mas- 
ter equation for quantum Brownian motion enabled us 
to obtain analytical expressions for the asymptotic en- 
tanglement and for the boundaries between the phases. 
In the present paper we complete and generalize the ideas 
presented in [24 . Here, we consider two different models 
for the coupling between the system and the environ- 
ment: We analyze not only the case where the coupling 
is bilinear both in the position of the system and the envi- 
ronment but also we solve the case where the coupling be- 
tween the system and the environment is symmetric be- 
tween position and momentum (this is technically equiv- 
alent to the RWA). We analyze in detail how the phase 
diagram changes depending on the coupling to the envi- 
ronment as well as on the environmental spectral density. 
For both models we study the entanglement for ohmic, 
sub-ohmic and super-ohmic spectral densities. Finally, 
we study the entanglement between non-resonant oscil- 
lators where a new master equation is derived. In such 
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case, we show that although non-resonant effect tend to 
eliminate entanglement, it is possible to have resilient 
entanglement at sufficiently low temperatures. 

The paper is organized as follows. In Section |IT] we 
review the basic technical tool we will use in our anal- 
ysis: the master equation. In Section |f 1 1 1 we show how 
to use this equation to analyze the evolution of the en- 
tanglement for general Gaussian states. In Section |IV| 
we present a detailed analysis of all qualitatively differ- 
ent evolutions (dynamical phases) of entanglement. In 
Section |V] we study the evolution of entanglement for 
non-resonant oscillators. In Section IVII we summarize 
and conclude. 



II. TWO EXACTLY SOLVABLE MODELS 

We will study the evolution of the entanglement be- 
tween two harmonic oscillators with coordinates xi and 
X2 (they constitute our system) which are coupled with a 
bosonic environment. We will analyze two different mod- 
els: First we will assume that the coupling between the 
system and the environment is bilinear in their position 
[221 123 HE]- Then, we will analyze the case where 
the coupling is symmetric in position and momentum. 
In both cases we will use an exact master equation to 
describe the evolution of the reduced density matrix of 
the system. In what follows we will briefly describe the 
two models and their solution. 



A. Quantum Brownian motion with position 
coupling 

The total Hamiltonian for the universe formed by the 
system and the environment is H = Hg + H^nt + -f^enu 
where 
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It is convenient to use coordinates x± = {xi ± X2)l\f^ 
since x+ couples to the environment. The Hamiltonian 
Hs is 

(p+ +P-) , m 2 2 , 2 2^, 
Am 1 

where the frequencies of the x± oscillators are — (wi + 
a;|)/2 ± C12 and the coupling constant between them is 

= — w|)/2. Below we will analytically solve a 

special but very important case: We consider the two 
oscillators to be resonant, i.e. we take Wi = L02 (in this 
case the x± oscillators are decoupled, as — 0). 



This model (known as Quantum Brownian Motion) 
can be exactly solved [28]. Thus only two parameters 
are necessary to characterize the effect of the environ- 
ment on the system. The first one is the initial state 
of the environment (assumed to be thermal, with initial 
temperature T). The second one is the spectral density 
of the environment, which is a function of the frequency 
defined as J(uS) = ~ Wn)/2'mnWn- One can 

show that the reduced density matrix p, obtained from 
the state of the universe by tracing out the environmen- 
tal oscillators, obeys an exact master equation which is 
written as [251125] : 

p = -i[HR,p]-i'y{t)[x+,{p+,p}]- 

~ D{t)[x+,[x+,p]]-f{t)[x+,[p+,p]]. (3) 

Here, the renormalized Hamiltonian is 

Hn = Hs+'^5u;'{t)xl. (4) 

The coefficients 6u!'^{t), j{t), D{t) and f{t) depend on the 
spectral density of the environment {D{t) and f{t) also 
depend on the initial temperature T). The explicit form 
of these coefficients is rather cumbersome and was stud- 
ied in detail elsewhere [251 [30] • Some results on the be- 
havior of the coefficients for typical environmental spec- 
tral densities will be described below. In particular we 
will consider the family of spectral densities of the form: 

J(c^) = ^m7oCc>(|)" '9{A-co), (5) 

where A is the cutoff frequency and 70 is a coupling con- 
stant. Depending on the value of n, the spectral densities 
are known as: ohmic (n = 1), sub-ohmic (n < 1) and 
super-ohmic (n > 1). 

To study analytically the long-time regime, we just 
need to assume (as it is the case for realistic environ- 
ments) that the coefficients of the master equation ap- 
proach asymptotic values after a temperature-dependent 
time. The time dependent frequencies fif 2 (*) ='^12 + 
6LLp'{t)/2 approach cutoff independent values only if the 
bare frequencies 0^1.2 have an appropriate dependence on 
the cutoff. The coupling constant ci_2 must also be renor- 
malized in the same way so that the time dependent cou- 
pling Ci2{t) = C12 + 5uP'{t)/2 approaches a finite cutoff 
independent value. The behavior of the diffusion coef- 
ficients D{t) and f{t) is more complicated and depend 
on the initial temperature. A word on notation: upper 
case letters will be used to denote renormalized quanti- 
ties. The time label will be omitted when referring to 
the asymptotic value of the corresponding function (i.e., 
f2i^2 denotes the asymptotic value of the renormalized 
frequency of the oscillators, etc). 

The master equation is a powerfuU tool to understand 
the qualitative behavior of the system. For this purpose, 
it is convenient to use it to obtain simple evolution equa- 
tions for the second moments of x± and p±. Thus, it is 
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simple to show that the second moments of x+ and 
satisfy the following equations: 
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Where ^l{t) is the renormalized frequency of the x+ os- 
cillator. In turn, the evolution equations for the second 
moments of a;_ and p_ are simply the ones of a free os- 
cillator (i.e., can be obtained from the above ones by 
considering vanishing values for all the coefficients of the 
master equation). 

From the above equations the interpretation of the co- 
efficients appearing in the master equation is transpar- 
ent: 7(t) is responsible for relaxation since it induces the 
decay of energy, D(t) is a normal diffusive term which 
induces heating increasing the momentum dispersion. In 
turn, f{t) the so-called anomalous diffusion coefficient is 
responsible for the squeezing of the asymptotic state (see 
below) or, in other words, of a violation of the equipar- 
tition principle: Thus, in the stationary state (which is 
reached only if the environment is such that the coeffi- 
cients approach constant asymptotic values) eq. (|7| im- 
plies that the expectation value of kinetic and potential 
energy differ by a factor which is proportional to fit). 
The role of this term, will be very important in our anal- 
ysis below. 

Our analysis will be based on the use of the above 
equations to study the long time regime for cases where 
the environment is such that the coefficients of the master 
equation do approach a constant asymptotic value. Thus, 
it will be useful to write down the explicit asymptotic 
values of the dispersions A^x^ — (a;^) and A^p+ = (p^). 
From the above equations we find that 
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and {{x+,p+}) = 0. It is worth noticing that depending 
on the sign of the asymptotic value of / the nature of 
the relation between the variances, or squeezing, may 
change quite dramatically. The sign of / indicates what 
observable is being effectively localized. In fact, if the 
coefficient / is positive the asymptotic state is localized 
in position (i.e., the equilibrium state is squeezed along 
position) , which is a feature of low temperatures. 



B. Quantum Brownian motion with coupling 
symmetric in position and momentum 

We will also consider another exactly solvable model 
which is very similar to the above one. The only differ- 
ence is that the system and the environment are coupled 



through different observables. The interaction Hamilto- 
nian between the two resonant oscillators and the envi- 
ronment is 
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In the case c„ — Cn the total interaction can be rewritten 
in terms of creation and anihillation operators of the x+ 
oscillator (denoted a and a^) and the environmental ones 
(denoted as &„ and blJ. Thus, 
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(afojj -I- a^n)- 



(10) 



This is the same type of interaction that one obtains 
by making the so-called rotating wave approximation 
(RWA) for the model with Hamiltonian It is worth 
pointing out that we will discuss this as a separate model 
with its own exact solution (and not necessarily as an 
approximation to the previous one). Here, the system 
interacts with the environment both through position 
and momentum. As in the previous case, interactions 
within the system are induced through the environment. 
Such interactions generate a renormalization of the sys- 
tem's parameters. To be able to properly renormalize 
all the parameters in the Hamiltonian of the system we 
should include the most general type of interactions in 
such Hamiltonian. In this case, it includes not only cou- 
pling between the oscillators coordinates but also mo- 
mentum coupling. Thus, the Hamiltonian of the system 
is 



Hs = 
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In the resonant case we are considering here, this Hamil- 
tonian is simply written in terms of coordinates x± as 
the sum of two decoupled oscillators with frequencies 
ujj. = ± ci2/w^)(l ± £12/0;^) and masses m± — 

m/(l±ci2V). 

In this case, it is possible to obtain an exact master 
equation for the reduced density matrix p. For the zero 
temperature case the exact master equation was obtained 
by An et al [31 . Their result can be generalized to finite 
temperature (details of the derivation will be presented 
elsewhere) and reads: 

p = -i[IlR,p\-i^{t)(^[x+,{p+,p}]-[p+,{x+,p}]^ 

- D{t)([x+,[x+,p]]-\ ^-^[p+,[p+,p]] 



Here, the renormalized Hamiltonian Hf( is 
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The main features of this master equation are sim- 
ple to understand: Not surprisingly, this equation looks 
as the symmetrized version of ([S]). Thus, the damping 
coefficient j(t) appears multiplying a term that is sym- 
metric under canonical interchange of position and mo- 
mentum. This is also the case for the normal diffusive 
term (proportional to D(t)). The absence of the anoma- 
lous diffusion is precisely an expected consequence of the 
same symmetry since this term is anti-symmetric in ([s]) . 
Renormalization is also symmetric since this type of cou- 



pling induces renormalization not only on the oscillator 
frequency but also on its mass. In fact, renormalized fre- 
quencies and masses of each oscillator can be defined as: 
n^{t) = Lu{l + Sn^{t)/2uj^), M,{t) = m/{l + 5n'^{t)/2uj'^). 
In turn, renormalized coupling constants are: Ci2(t) = 
ci2 + <5fi2(i)/2, Ci2{t) = 5i2 + sn\t)/2. 

From the master equation we can obtain equations of 
motion for the second moments of the oscillator: 



{{x+,p+})-A^{t){xi) + --——D{t), (13) 



where M{t) = m/(l -I- {SCfit) + cxi)lu?) and VL(t) = 
Lo(\ -I- ((5f22(i) -|- c\2)IliP') are the mass and the frequency 
of the oscillator xj^. The role of each term is transparent: 
7(i) is an effective damping rate inducing decay towards 
the ground state while 13 (t) is a symmetrized diffusion 
constant inducing the spread of the state both in position 
and momentum. Assuming these coefficients approach 
constant asymptotic values we can easily derive the long 
time values of position and momentum dispersions to be 
given by 



(14) 



and = 0. 

Contrary to what happened in the non-symmetric case, 
governed by the master equation ([3|, the asymptotic 
state satisfies the equipartition principle since the expec- 
tation values of kinetic and potential energies are identi- 
cal. Analogously, as will be mentioned below, the asymp- 
totic state of the a;+ oscillator is not squeezed. 



III. EVOLUTION OF ENTANGLEMENT FOR 
GAUSSIAN STATES 

We will assume that the initial state of the system is 
Gaussian. As the complete evolution is linear, the Gaus- 
sian nature of the state will be preserved for all times. 
This enables us to analytically compute the entanglement 
between the two oscillators in the following way: Entan- 
glement for Gaussian states is entirely determined by the 
properties of the covariance matrix defined as 



where i,j = 1,...,4 and r = {xi,pi,X2,P2)- In fact, a 
good measure of entanglement for such states is the so- 
called logarithmic negativity Ej^ [32l [33] which can be 
computed as P^ I551[M) : 



Ej^ = max{0, - ln(2i^min)}, 



(16) 



where Vrnin is the smallest symplectic eigenvalue of the 
partially transposed covariance matrix. There are known 
expressions for Ej^ for particularly relevant Gaussian 
states which will be used as initial conditions in our study. 
For this reason it is useful to mention them here: For the 
two-mode squeezed state, obtained from the vacuum by 
acting with the operator exp(— r(a|a2 — 0102)), we have 
Ej^ = 2\r\. For this state the dispersions satisfy the min- 
imum uncertainty condition Sxj^^Sp^ — Sx-Sp- = 1/2. 
The squeezing factor determines the ratio between vari- 
ances since mflSx^ / Sp-^- = 6p-/{mflSx^) — exp(2r). As 
r — > cx) the state becomes localized in the p^ and a;_ 
variables approaching an ideal EPR state [35j . 

Now we consider a general initial gaussian state of the 
two oscillators. From the appropriate master equation 
([3]) or ( 1 1 ) we showed how to obtain equations for the co- 



V.,{t) = {{r,,r,})/2~{n){r,), 



(15) 



variances. These equations split into two blocks of 2 x 2. 
The evolution of the first block formed with the second 
moments of .t_ and p_ corresponds to a free oscillator 
with frequency which can always be expressed in 
terms of two dispersions Sx^ and 6p^. The evolution 
equations for the second block, formed with the second 
moments of x^ and p^, were discussed above and yield 
equilibrium values Aa;_|_ and Ap+. It can be also easily 
proved that the off-diagonal block, containing the corre- 
lations between the oscillators (x+,a;_), vanishes in the 
asymptotic regime. These simple observations are al- 
most all we need to fully analyze the evolution of the 
entanglement between initial Gaussian states. Thus, us- 
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ing the diagonal block form of the covariance matrix in 
the (a;-|-,a;_) bases (and changing basis to obtain covari- 
ances of the original xi^2 oscillators) it is simple to find 
the smallest symplectic eigenvalue of such matrix and 
compute the logarithmic negativity. The result is: 



where the function E{t) is defined as 

E{t) = EM + ^EMG{t). 



(17) 



(18) 



Here G{t) is an oscillatory function with period ■k/uj^ 
that takes values in the interval {—1,+!}. Its explicit 
form will be given below. The mean value Ej^ and the 
amplitude lS.Ej^ that characterize the oscillations of E{t) 
are simply written as 



E^f = max{|r|, |rc„t|} 
l^Ef^ = min{|r|, \rc„t\}- 



(19) 
(20) 



In the above equations r is the initial squeezing factor 
defined as 



r = - In 
2 



5x- 



(21) 



and Tcrit is related to the squeezing factor of the equilib- 
rium state for the x-L-oscillator 



rcrit = 2 



Finally, Scrit is defined as 



Aa 



Ap+ 



(22) 



Scrit = - ln[4:Ax+Ap+Sx-Sp-], 



(23) 



and turns out to be simply related with the entropy of 
the asymptotic state. Thus, the von Neumann entropy of 
the final state (S^) is Sy = f{<J+) + f{(J-) where f{a) — 
((T+i)ln(cr+i)-(CT-i)ln(cr- i), with a+ = Ax+Ap+ 
and cr_ — Sx-dp-. It is worth mentioning that in all the 
above formulae the dispersions Aa;+ and Ap_|_ are the 
asymptotic values of the dispersions along position and 
momentum that depend upon the temperature and the 
type of coupling to the environment (and that, for the 
models analyzed above are given by eqs. Q and (14)). 
For completeness we give the explicit formula for G{t) 
which is, indeed, not very illuminating: 



AE^G{t) = max{|r|, \r^„t\} + i^ln 



cosh[2(r — Tcrit)] cos^{uj^t) + cosh[2(r + rent)] sm'^{uj^t) 



\j2[smh^{2r) + sinh^(2rc,.it)) sin^(u;_t) cos2(tj_t) + sinh^[2(rcrit - r)] cos4(u;_t) + sinh^[2(rcrit + r)] sin^(a;_t) 



These simple results will enable us to draw the fol- 
lowing conclusions about the dynamics of entanglement 
for long times. 1) Evolution of Ej^ is fully character- 
ized by Tcrit and Scrit- 2) Only three qualitatively dif- 
ferent scenarios emerge. First, entanglement may per- 
sist for arbitrary long times. This phase, which we call 
"NSD" (for no-sudden death) , is realized when the initial 
state is such that Ej^/ — AEf^ > 0, which translates into 
Ikl — kcritll > Scrit- Then, there is a phase where entan- 
glement undergoes an infinite sequence of events of "sud- 
den death" and "sudden revival" [37] . This occurs if 
the initial state is such that \Ec\ < r < —Ec + 2\rcrit\, 
where the quantity E,. is defined as 



Ec 



S f 



(24) 



We denote this phase as "SDR" (for sudden death and 
revival). Finally, a third phase characterized by a fi- 
nal event of "sudden death" of entanglement may be re- 
alized if |r| < —Ec- This phase is simply denoted as 
"SD" (for sudden death). In what follows we will analyze 
these phases for different spectral densities and coupling 
strength between the oscillators and the environment. 



Some further physical insight about the origin of the 
entanglement can be obtain by rewriting the eq. ( 18 1 as 



E{t) = \rcr^t\ - Scr^t + \r\G{t), if |r| < \rcr^t\, 
E{t) ^ \r\ - Scr.t + \rcrit\G{t), if \r\ > \rcr^t\- 



In this way it is clear that for initial values |r| < \rcrit \ the 
environment supplies the resource to generate entangle- 
ment. In particular, when jrcr-itl — Scrit ^ 2|r|, then the 
entanglement in the final state is larger than the quan- 
tum resource (squeezing) available in the initial state. In 
other cases the environment does not act as the supplier 
but simply degrades the quantum resource which is al- 
ready present in the initial state (either in the form of 
squeezing or entanglement). Below, we will analyze this 
further by using a very convenient tool: a phase dia- 
gram where the fate of entanglement can be graphically 
depicted for all initial states. 
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IV. EVOLUTION OF ENTANGLEMENT 
A. Phase diagrams for entanglement dynamics 

Here we will introduce a convenient tool to study the 
different dynamical phases of entanglement. In fact, de- 
pending on the properties of the environment (initial 
temperature, damping rate, etc) a given initial state (pa- 
rameterized by the squeezing r and by the product of ini- 
tial dispersions 5x^5pJ) will belong to one of the three 
phases: SD, NSD or SDR. For fixed values of 7 and 
Sx^dp^ we can always draw a phase diagram like the one 
displayed in Fig. [T] To obtain it we need to analyze the 
temperature dependence of the asymptotic dispersions to 
obtain both \rcrit\ and Scrit as a function of the temper- 
ature. In the phase diagram, the areas corresponding to 
each of the three phases are displayed. As a reference, 
we also include two curves that show the temperature 
dependence of Scrit and \rcrit\ (dashed and dotted lines 
respectively). The actual diagram shown in Fig. T] corre- 
sponds to a particular case: an environment with ohmic 
spectral density coupled to the system through position 
with C12 = (then, = H. and r7i_ = m). We also 
assumed a pure initial state with 5x_dp^ = 1/2. In the 
following sections we will see that other spectral densities 
will give rise to slightly different features in the phase di- 
agram but its topology will remain unaffected. Changes 
in the initial state (i.e., considering mixed states with 
Sx-dp- > 1/2) can also be simply understood and will 
be discussed below. 

The phase diagram describes all dynamical informa- 
tion of the asymptotic evolution of entanglement for the 
case of position coupling to the environment (see below 
for symmetric coupling) . Some important features of the 
phase diagram are worth mentioning. In particular, we 
would like to focus first on the NSD phase present at low 
temperatures. Its origin is purely non-Markovian and 
non-perturbative. Its area shrinks as the damping rate 
decreases. The states in this phase are the ones for which 
the final entanglement may be larger than the squeezing 
invested in the initial state. For such states the entangle- 
ment mostly comes from the squeezing available in the 
environment. This is particularly clear for the case of 
coherent states, that can become entangled below the 
critical temperature Tq (see below). 

To understand the nature of entanglement in this re- 
gion of the phase diagram it is useful to focus first on 
the properties of such diagram along with of its axis. 
The zero temperature line (i.e., the horizontal axis) con- 
tains states in the NSD phase for small and large squeez- 
ings. Thus, the NSD phase is realized at zero temper- 
ature when the initial squeezing r is either \r\ > r-i or 
Irl < ri, where 
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2TnVLl\x\{T = ^) 
2Apl{T = 0) 



(25) 
(26) 



see Fig. [T] For the range of squeezings between ri and 
r2 (the region centered about |rcrit|) the states belong to 
the SDR phase. This implies that pure initial states (at 
T — 0) will never experience a sudden death. They will 
never belong to the SD phase. 

It is interesting to notice that for T = the asymp- 
totic state of the x^ oscillator is squeezed in position (i.e., 
m,nAx+{T = 0) < Ap+(T = 0)) and that it also has non- 
vanishing entropy (i.e., Ax+{T = 0)Ap+(T = 0) > 1/2). 
In the following section we will present analytic expres- 
sions for Aa;_|_ and Ap_|_ in the case of the ohmic environ- 
ment. Here it is sufficient to mention that the squeezing 
in position is a consequence of the fact that, being the 
interaction with the environment through the position 
observable, the asymptotic state tends to localize more 
along position than along momentum. The squeezing ri 
(r2) is precisely the ratio between the asymptotic posi- 
tion (momentum) dispersion and the one corresponding 
to the vacuum: a non-vanishing value of ri means that 
the asymptotic state of the oscillator x+ has a dispersion 
along position that is smaller than the vacuum. There- 
fore, the states that belong to the low temperature NSD 
island are the ones for which the state of the oscillator 
x^ is narrower in position than the vacuum dispersion 
which is given by l/2mfi. As the temperature of the 
environment increases, the asymptotic value of the po- 
sition dispersion Aa;+ also grows. Therefore, the NSD 
phase shrinks and completely disappears above the crit- 
ical temperature Tq, which is precisely the one for which 
the position dispersion becomes identical to the one cor- 
responding to the vacuum, i.e. 



To such that Aa;+(T = Tq) 



(27) 



The vertical axis of the phase diagram it is also worth 
analyzing since it describes the fate of coherent states 
(for which r = 0). As mentioned above, for temperatures 
lower than Tq such states end up entangled due to the 
interaction with the environment. But for temperatures 
larger than Tq such states always experience an event of 
sudden death of entanglement (the states belong to the 
SD phase). 

The high temperature region of the diagram is rather 
different than the low temperature one. Thus, for high 
temperatures we have E,. < (which implies that coher- 
ent states do not get entangled) and also Vcrit Scrit 
(which implies that the region covered by the SDR phase 
becomes relatively narrower). Hence, initial states with 
large squeezing (|r| > ln(2Aa;_|_Ap_|_)/2 = Scrit) retain 
some of their entanglement while those with squeezing 
factors smaller than the critical value Scrit suffer from 
sudden death. However, our analysis shows that the 
boundary between SD and NSD phases is rather subtle: 
for any finite temperature the two phases are separated 
by a very narrow portion of SDR phase (in this phase 
there are oscillations of the entanglement whose ampli- 
tude, |rcrit|, depends on the temperature in a way which 
is different for different spectral densities, as will be dis- 
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cussed below). In any case, these oscillations are, indeed, 
yet another interesting non-Markovian effect identified by 
our analysis. 

A final comment on the phase diagram: The NSD 
phase is characterized by a non- vanishing asymptotic en- 
tanglement that can be quantified in a straightforward 
way from the phase diagram itself. The average value 
of the logarithmic negativity is simply the distance to 
the dashed line (which signals the midpoint of the SDR 
phase) or just the distance between the dashed and dot- 
ted lines for \r\ < \rcrit\- 



states (with r = 0) belong to the SD phase at zero tem- 
perature. For a given value of the temperature T, the 
asymptotic state of has some degree of mixing. The 
states that have enough squeezing r to support some 
entanglement for such degree of mixing are denoted as 
GLEMS (Gaussian least-entangled mixed states) |34|- As 
a final comment we should point out that for an initial 
two-mode squeezed state, the curve defining the bound- 
ary of the NSD region (given by eq. (28 1 ) coincides with 



the one obtained in |19j , where entanglement was studied 
under a Markovian rotating-wave approximation. 




FIG. 1; Phase diagram for ohmic environment {ft = I, m = 1, 
70 = 0.1, A = 20, Ci2 = 0, 5x-Sp- = 1/2). The sudden 
death (SD), no-sudden death (NSD) and sudden death and 
revival (SDR) phases describe the three different qualitative 
long time behaviors for the entanglement between two oscilla- 
tors interacting with the same environment. The SDR phase 
is centered about the dashed line Scrit and has a width given 
by the dotted line |rcrit|. This is the case for temperatures 
above To, the one for which Scrit = Ij'critl. Below this tem- 
perature the role of Scrit and \rcrit\ are interchanged. SDR 
separates the SD and NSD phases. The low temperature NSD 
island is due to non-Markovian and non-perturbative effects. 
Ej^ in the NSD phase is the distance to the dashed line for 
\r\ > \rcrit\, and the distance between the dashed and dotted 
lines for \r\ < \rcrit\- 

There is an important qualitative difference between 
the cases of position coupling and symmetric coupling. In 
the later case, the symmetry implies that the asymptotic 
state of the x-^ oscillator is not squeezed as Vcrit = 0. 
The asymptotic entanglement is Ej^{t) = max{0, i?(t)}, 
where: 



- ln[4A2;_|_Ap-|_(5a::_(5p_]. 



(28) 



This implies that for symmetric coupling the phase 
diagram is simpler, as shown in Fig. [2] In this case 
there are only two phases (SDR does not exist). The 
NSD phase is characterized by the condition \r\ > 
ln{4Ax+Ap+Sx-Sp-)/2. In this case, the entanglement 
achieved in the asymptotic regime is the difference be- 
tween \r\ and the curve that limits the two phases. Con- 
trary to what happens for position coupling, initial pure 




0.5 1 1.5 2 
Irl 

FIG. 2: Phase diagram for ohmic environment with symmet- 
ric coupling (n = 1, M = 1, Ci2 = 0, SxSp- = 1/2). The 
phase diagram is qualitative different from the correspond- 
ing to the position coupling. The SD and NSD phases are 
present. Vcrit ~ for every temperature and the asymptotic 
entanglement is allways constant. Here E^^ is the distance 
from \r\ to the line that limits the two phases. 



Evolution in the different phases: Analytic and 
numerical results. 



Here we will analyze the above results contrasting the 
analytic predictions with the results of an exact numer- 
ical solution of the problem. Numerical solution is in- 
deed exact since in the case of a discrete environment 
(formed by N oscillators). It is obtained by using the 
linearity of the problem to exactly evolve the complete 
covariance matrix and to obtain the full quantum state. 
Once this is done one can directly compute the logarith- 
mic negativity (see [38l El] for another application of 
the same method) . Whenever possible (position coupHng 
with ohmic spectral density) we compared this evolution 
with analytic expressions for the exact reduced evolu- 
tion operator, finding complete agreement between both 
methods. 



Position CO 



a) Ohmic spectral density: Now, we will focus on the 
ohmic environment eq. ([3| {n = 1), where the high fre- 
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quency cutoff A defines a characteristic timescale 
over which the coefficients j{t) and 6uj'^(t) vary. For 
times t » A^^ these two coefficients settle into asymp- 
totic values: ^{t) — > 7 = 270 and Suj'^{t) ~4A'^/tt. It 
is worth mentioning a technical point related with the 
renormalization that seems to have caused some confu- 
sion in the literature. The interaction with a common 
environment induces a coupling between the oscillators. 
Thus, even if we consider a vanishing "bare" coupling 
(i.e., C12 = 0) the asymptotic value of the coupling will 
be non-zero and given by C12 — Suj'^/2. It is natural to 
define renormalized parameters of the oscillators as the 
ones characterizing the long time limit. Thus, for the 
renormalized coupling to be C12 = we must consider a 
bare coupling C12 = — 5w^/2 in the original Hamiltonian. 
This simply says that the coupling constant between the 
oscillators must be renormalized in the same way as their 
natural frequency (with the same counterterm) . If one 
does not do this (and assume, for example, that the bare 
coupling vanishes) one would observe high frequency os- 
cillations at long times (with a frequency which is set by 
the cutoff A) . On the contrary, by adding the appropri- 
ate counterterms to the bare Hamiltonian one obtains a 
A-independent long time limit. In such case, we have 
fii = fi2 = = i-^-- 

Predictions discussed in the previous Sections can be 
verified by an exact numerical solution to the problem. 
For our numerics we considered parameters 70 = 0.1, 

= 1, A = 20, m — 1, C12 = (extension to the 
case where the natural oscillators interact can be easily 
done) . We considered separable squeezed states for which 
mfifei 2/<5pi,2 — exp(2r) as well as two-mode squeezed 
states for which mfl6x+ /^P+ — Sp^/(mflSx-.) — exp(2r) 
(in both cases Sx^Sp^ = 1 /2). In Fig. |3] we show the en- 
tanglement dynamics in an environment at zero tempera- 
ture. We clearly see that the final entanglement achieved 
by different initial states only depends upon the squeez- 
ing factor r. Initial entangled states reduce their degree 
of entanglement while initial separable states do get en- 
tangled through the interaction with a common environ- 
ment. Evolution of separable states with positive and 
negative squeezing is compared in Figs, [s] (a) and (b). 
In the first case entanglement grows much faster. This is 
due to the fact that the initial state has a wider spread in 
the position observable, which is the one appearing in the 
interaction Hamiltonian. In the asymptotic regime, as it 
is predicted, the dynamics is the same. They oscillate 
with the same frequency around the same mean value 
and with the same amplitude, but as it was expected, 
they have a phase shift of 7r/2. 

The existence of events of sudden death and revival can 
also be seen from the numerical solution and are shown 
in Fig. |4](our numerical results show full agreement with 
the analysis presented above concerning the nature of the 
SDR phase). In the same Figure, we also show the evo- 
lution belonging to the NSD phase. They correspond to 
a squeezing such that \r\ < \rcrit\- In such case the am- 
plitude of oscillations in the asymptotic regime is equal 



to |r| and the mean value is Ec- An example of the SD 
phase appears in Fig. [4] along with another example of 
the NSD phase for a non-zero temperature. It is also 
noticeable that the amplitude of the oscillations almost 
vanishes in the high temperature limit. 
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FIG. 3: Logarithmic negativity for resonant oscillators in the 
same environment, (a) For T = the NSD phase appears 
both for large and small squeezing. Initially separable states, 
squeezed or coherent can get entangled. The asymptotic be- 
havior only depends on r. The amplitude of oscillations van- 
ishes when r — > 0. (&) Initial states with negative squeezing, 
have the same asymptotic behavior with a dephasing of 7r/2. 
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FIG. 4: (a) The SDR phase appears for intermediate values of 
squeezing at zero temperature. For |r| < \rcrit\ (dashed line) 
E]^ = Ec and the amplitude of oscillations is equal to \r\. 
[b) T/Q, = 10, the SD phase appears for small |r| and NSD 
phase for large squeezings, oscillations in the steady state are 
attenuated as the temperature increases. 

We can obtain an analytical expression for the parame- 
ters needed to analyze the entanglement dynamics. Thus, 
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using the exact expressions obtained in f3IP we find that 
at zero temperature 



n = E,{T = 0) = -In 



r2 



In 



2 - Aj^/n'^ 

v/1 - ly^' 



2 arccos(7/ri) 
arccos(7/ri) 



4 7 
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Scrit — 4 



+ 27/r!v/r^7W 



hi[A/r2] 
arccos(7/ri) 



4 1 ~ 27V»^ 
^ 1 -72/f^2 
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arccos(7/ri) 



These formulae have a simpler form in the weak cou- 
pling limit where 
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In this case, the asymptotic coefficients of the master 
equation up to second order in 7 are given by: 



D 



f 



2mj'^ 



'2ln 



rA 

n 



27 



In 



(29) 



A technical comment is in order here: To estimate the 
asymptotic behavior using an expansion in powers of the 
coupling constant, we need the coefficient D to one order 
higher than /. This fact was already noticed in [3D] and 
is evident from the fact that critical squeezing is given 

by 



1 



'^crit 



■In 



1 



2to7/ 



D 



(30) 



Estimates for the critical temperature Tq (the temper- 
ature for which the position dispersion becomes identical 
to the vacuum one) can be obtained as follows: position 
dispersion at low temperatures is: 



T 



1 



zim 



H 



'1- 



■r 



27rT 



where the function H{z) is the Harmonic Number. Ex- 
panding this for low temperatures (T/fl <C 1), 



arccos(7/ri) 



15 



(l-2(7/^^)^)(^) 



(31) 



we can obtain an approximate expression for 
To that accurately reproduces our results for 
arctan (^1-7/^2)2/(7/^2)) < n/2. 

On the other hand, for high temperatures we can use 
the appropriate approximations to obtain: 
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27 
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A + n 



Scrit ~ 2 
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27 I A + n 



■In 



A-n 



In this regime Vcrit approaches a temperature- 
independent value that decreases with the high frequency 
cutoff and increases with the coupling constant 7. As 
a consequence, the asymptotic entanglement is approx- 
imately constant. The behavior of Scrit is simpler: as 
expected it behaves as the entropy, growing as ln(T) for 
high temperatures. The narrow passage between the SD 
and the NSD phases closes as 1/A and moves to larger 
and larger values of squeezings. For completeness we in- 
clude the diffusion coefficients in the high temperature 
regime up to first order in 7. They are 



D 



f 



2m'jT, 



27 



In 



A + n 
A-n 



T. 



(32) 



h) Sub-ohmic spectral density: Here we will analyze 
the behavior of entanglement in an environment with 
a sub-ohmic spectral density as ([s]) with n = 1/2. In 
this case, the oscillators of the infrared and the resonant 
bands are coupled more strongly to the system (since 
f2 ^ A) and the environment induces more dissipation. 
As a consequence, the equilibrium state of the oscillator 
is noticeably more squeezed along position than the 
one corresponding to the ohmic case |40j . Therefore, con- 
sidering eq. (20 1, we expect a larger value for \rcrit\ at 



zero temperature which, in turn, would imply that the 
oscillations of the entanglement in the steady state will 
have larger amplitude. In addition, the entanglement for 
initial coherent states will be larger as well as the critical 
temperature Tq. Below we will show only the numerical 
results, since there are no available analytic expression 
for the coefficients of the master equation. We used the 
same parameters as in the previous subsection, noticing 
that in this case Suj'^{t) — 8(27o)/7rA. 

In Fig. [5] we show the dynamics of entanglement for 
two resonant oscillators immersed in a bath at zero tem- 
perature. There we can appreciate oscillations of entan- 
glement in the steady state with larger amplitude than 
in the ohmic case but with the same frequency. For ini- 
tial coherent states, the entanglement achieved is greater 
than in the ohmic case. As mentioned above, this is a 
consequence of the coupling between the system and the 
resonant bands of the environment that produce a sub- 
stantial squeezing in the steady state. We can also notice 
that the system approaches equilibrium earlier than in 
the ohmic case due to the fact that dissipation is stronger 
than in the ohmic case. 
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FIG. 5: Logarithmic negativity for resonant oscillators in a 
sub-ohmic environment, (a) For T = the NSD phase ap- 
pears both for large and small squeezing. The amplitude of 
oscillations are higher than in the ohmic case. And as a con- 
sequence the entanglement achieved for initial coherent states 
is also higher. (6) For T/fi = 10 appreciable oscillations are 
present in the asymptotic regime. 



The asymptotic features of the entanglement observed 
in the figures for specific cases, can be summarized in 
the phase diagram shown in Fig. [6] The shape of this 
diagram is essentially the same as the preceding case. In 
the low temperature regime we find again an NSD island 
with an area that is larger than the one corresponding 
to the ohmic case. Also the value of |rcrit| at zero tem- 
perature is larger than for the ohmic environment, and 
it decreases with the temperature. As a consequence, 
we also observe a SDR region at high temperature whose 
width is given by |rcrit|. In this case, the oscillations that 
appear in the ohmic spectral density are enhanced. 




FIG. 6: Phase diagram for sub-ohmic environment. The SD, 
NSD and SDR phases that describe the three different qual- 
itative long time behaviors for the entanglement are present. 
The low temperature NSD island is bigger than the one cor- 
responding to the ohmic environment. 



Even though there are no analytic expressions for this 
environment, we can obtain approximate formulae in the 
weak coupling and high cutoff regime. We can use them 
to get some qualitative understanding of the expected 
behavior but they are useles to draw quantitatively con- 
clusions (since we can only obtain all coefficients of the 
master equation up to first order in 70 but, as indicated 
above, to estimate quantities such as Vcrit and Scrit we 
would need to have the asymptotic form of the coefficient 
D to one order higher). Thus, to lowest order we find 
7(t) — + ^suh — 27o\/A/ri, which grows with the cutoff as 
expected. At zero temperature, the diffusion coefficients 
are 



D sa mjsub^ 

f ^ 1: 



In 



A + n 
A-n 



(33) 
(34) 



Again, the two coefficients grow with the coupling and 
the cutoff frequency. The anomalous diffusion coefficient 
/ is larger than the one corresponding to the ohmic case. 
This produces a stronger localization of the asymptotic 
state in the position observable. We can also obtain the 
high temperature expressions in the weak coupling limit: 



D « 2m^subT, 
T 



f 



-27. 



(35) 
(36) 



In this case all the coefficients are proportional to the 
temperature. The fact that the asymptotic state is 
squeezed at high temperatures is a consequence of the 
fact that the coefficient / approaches significantly higher 
values than the ones corresponding to the ohmic case. 

c) Super-ohmic spectral density: A super-ohmic envi- 
ronment has a spectral density characterized by a higher 
population of high frequency bands. A typical example is 
given by equation (|5| with n = 3. Super-ohmic environ- 
ments are weakly dissipative. In fact, in this case the dis- 
sipation coefficient approaches an asymptotic value given 
by lit) — > Jsup — 27o(ri/A)^. The frequency shift is 
6u!'^{t) — > — 4(27o)A/37r. Thus, dissipation strictly van- 
ishes in the infinite cutoff limit. In such case the os- 
cillator a;_(. does not reach equilibrium, a fact that was 
also noticed in |41j and is related with the phenomenon 
of recoherence that could be induced by this type of en- 
vironment (i.e., decoherence is reversible in this case). 
In this limit we cannot apply the analysis presented in 
the previous Sections, which requires the oscillator x+ 
to approach equilibrium. Thus, for a super-ohmic en- 
vironment we expect to observe an oscillatory behavior 
for the entanglement up to very long times. This is pre- 
cisely what is observed in Fig. [7] where we show the 
results of the numerical solution for two different initial 
states. Oscillations of entanglement persist for low and 
high temperatures. The amplitude of the oscillations de- 
creases very slowly with time simply due to the fact that 
we consider a finite value for the high frequency cutoff. 
For this reason the value of the dissipation coefficient is 
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not strictly zero but very small. This implies that the 
system would reach an equilibrium in the extremely long 
time limit (i.e. for times of the order of es- 
timate which is consistent with the numerically observed 
behavior) . It is worth noticing that the result we present 
here for the super-ohmic environment is not compatible 
with the ones reported in [22' (the super-ohmic results of 
that paper seem to be simply in error, a more detailed 
comparison with such results will be presented below). 
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FIG. 7: Logarithmic negativity for resonant oscillators in a 
super-ohmic environment with 70 = 0.15. (a) For T = we 
do not observe that the entanglement achieves equilibrium. 
Oscillations are present for long times, (b) At T/fi = 10 the 
entanglement oscillates with smaller amplitude. 

For zero temperature we can also obtain the asymp- 
totic behavior of the diffusive coefficients in the weak 
coupling limit. These coefficients behave as f ~ 27o/7r-|- 
7s„pln[(A2 - f)2)/f]2j^^ ^ 27o/7r and D « m^jsup- 
Thus, in this case the anomalous diffusion / is propor- 
tional to the coupling constant and becomes independent 
of the cutoff. It takes the smallest value, comparing the 
three spectral densities that we considered, which is a sig- 
nature of the weak coupling between the system and the 
resonant band of the environment. On the other hand, D 
vanishes in the infinite cutoff limit (as mentioned above, 
Jsup vanishes as well). In the high temperature regime 
we have / « 27oT/7rA and D « 2mT'^sup- Here the 
small value of / produces a squeezing of the asymptotic 
state which is smaller than the one achieved for ohmic 
and sub-ohmic environments. 



2. Coupling symmetric in position and momentum 

Here we will consider the case where the coupling to 
the environment is symmetric in position and momen- 
tum. This model at zero temperature was studied pre- 
viously in f22J with a two-mode squeezed state as the 



initial condition. Here we extend these results by consid- 
ering arbitrary initial Gaussian states and arbitrary tem- 
peratures (we also take the opportunity to correct some 
erroneous results reported in [52^). The main conclu- 
sion concerning entanglement dynamics was announced 
before: only two phases (NSD and SD) exist. This con- 
clusion follows from the fact that the master equation is 
symmetric under canonical interchange between position 
and momentum. It is independent of the precise form of 
the asymptotic values of the coefficients appearing in the 
equation as long as equilibrium exist (which is not the 
case for the super-ohmic environment). 

We confirm this by a detail study of the numerical solu- 
tion using the same parameters of the previous sections. 
Here we also considered C12 — C12 — (then, Q, = lj_ 
and M — m__). In Fig. |8]we show the dynamics of en- 
tanglement for ohmic and sub-ohmic spectral densities. 
Our result show that at zero temperature entanglement 
is reduced to exactly half of its initial value [22J. This 
is a prediction of eq. (19 1 which is valid both for ohmic 



and sub-ohmic spectral densities, since in all those cases 
the asymptotic state of the oscillator x+ is the ground 
state. In fact, the form of the master equation at zero 
temperature ensures the stability of the ground state. 
As the asymptotic state is pure then this process can 
be thought of as a way to create pure gaussian entan- 
gled states from initial separable ones. Another obvious 
consequence of the symmetric coupling is that the behav- 
ior of initial states with negative or positive squeezing is 
identical. In Fig. [8] we show an example of the behavior 
of entanglement at temperature different from zero. As 
we discussed above, the steady state has non-zero en- 
tanglement. These results are summarized in the simple 
phase diagram of Fig. [2] which is essentially the same 
both for ohmic and sub-ohmic environments. 

For a super-ohmic environment the dissipative coeffi- 
cient 7 scales as 1/A^. In Fig. [9] we show that entan- 
glement oscillates for very long times, a simple conse- 
quence of the vanishingly small value of the dissipative 
coefficient. This contradicts the results obtained in [22] 
where it was shown that for a super-ohmic environment 
entanglement achieves equilibrium before the ohmic and 
sub-ohmic cases. From our previous analysis, based on 
the use of the master equation, we can simply conclude 
that the results of [25] do not seem to be reliable. On the 
contrary, our numerical results support the simple con- 
clusion obtained analytically by using the master equa- 
tion, which is local in time. Thus, entanglement oscil- 
lates slowly decaying with a rate that is roughly given by 
7 (which goes to zero in the infinite cutoff limit). 

It is simple to obtain analytic estimates for the asymp- 
totic values of the time dependent coefficients using a 
perturbative approximation. In that case, the dissipa- 
tion constant is 7 ^ 47o(ri/A)"^^ = 2J(ri)7r/Om. This 
shows that sub-ohmic environment induces stronger dis- 
sipation than an ohmic one. Also, for supra-ohmic en- 
vironments (n > 1) there is not equilibrium in the in- 
finite cutoff (Markovian) limit. In the same way, we 
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FIG. 8: Entanglement dynamics for resonant oscillators in 
an environment with symmetric coupling, (a) Environment 
at zero temperature, the asymptotic entanglement depends 
on the squeezing r and is constant. Ohmic and sub-ohmic 
environments arrive at the same equilibrium entangled state, 
(b) Environment at T/Q — 10, the final entanglement depend 
upon the initial squeezing for both spectral densities. 



can obtain the diffusion coefficient which is given by 
D = 2 J(ri)7rcoth(^). Using this expression the asymp- 
totic dispersions for the x+ oscillator are: 



9 9.9 .9 Mfl / ^ \ 
f2o2A„2 _ A„2 _ cothi 1- 



\2TJ' 



(37) 



These approximate expressions enable us to recover the 
results reported in jl9j . 




FIG. 9: Oscillators with symmetric coupling immersed in a 
super-ohmic environment. T = 0, the entanglement oscillates 
approaching approximately its initial value. The amplitude of 
oscillations decreases slowly since we are considering a finite 
cutoff. We also observe oscillations of entanglement for higher 
temperatures. 



V. NON-RESONANT OSCILLATORS 

The above properties are valid under a single impor- 
tant assumption: the two oscillators are resonant. If 
this is not the case the analysis becomes more compli- 
cated. The master equation is no longer valid since the 
x± modes are coupled. As X- is not isolated it also ap- 
proaches equilibrium. To analyze this we can obtain a 
new perturbative master equation (assuming the inter- 
action is through the position). It reads: 



-i[HR,p] 
.m 



ilit)[x+,{P+,P}] ~ D{t)[x+, [x+,p]] - fit)[x+, [p+,p]] 



- ^^S^+-it)[x+,U~,P}] -^il+-{t)[x+,{P-^P}] - D+-{t)[x+,[x^,p]] - f+^{t)[x+,[p-,p]]. 



(38) 



r 



As seen in the above equation, the x± oscillators inter- 



act with a coupling constant c_| ~ {loi 



\)I2 while 



a;+ is directly coupled to the environment. One of the 
terms coupling x± in ( 38 1 is a renormalization of the 
couphng constant. There is also a diffusive and a dis- 
sipative term. All the coefficients labeled with ± in- 
dices are proportional to the detuning A = (o^i — 0^2) 
(all of them vanish in the resonant limit). We can obtain 
asymptotic dispersions of the two oscillators, but the cor- 
responding formulae are rather cumbersome. The most 
important generic conclusions that we can draw from the 
above non-resonant master equation are the following. 1) 
As there is a final equilibrium state for both ± oscilla- 
tors, the final entanglement becomes independent of the 
initial state. 2) The approach to equilibrium proceeds 
with two different timescales (one of the decay rates is 



proportional to the detuning). 3) For sufficiently high 
temperatures the generic fate of the asymptotic regime 
is SD and the final state contains no entanglement. 4) 
However, for very low temperatures the final state can be 
entangled. The origin of the final entanglement resides 
again in the squeezing of the equilibrium state. If the 
± oscillators reach a final state with different squeezing 
(the squeezing of the a;+ oscillator is different from that 
of x_ mode) then the final state may be entangled. For 
very low temperatures this condition may be verified and 
entanglement may be present in the final state. In Fig. 
|10| we show how things change when we move away from 
the resonance condition if the environment is in a finite 
temperature state {T/Q. = 10). There we show the dy- 
namics for ohmic and sub-ohmic spectral densities. For 
early times, initial separable states become entangled due 
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It is interesting to notice that for any detuning it is 
possible to find temperatures below which the asymptotic 
state is entangled. The origin of this asymptotic entan- 
glement, as mentioned above, lies on the final squeezing 
of the x± modes. The dependence of the final entangle- 
ment with temperature is analyzed in Fig. |11| The curve 
is reminiscent of a phase transition with critical temper- 
ature depending upon the detuning. The existence of 
entanglement in the asymptotic state is not really a total 
surprise and is clearly related to recent findings of the ex- 
istence of entanglement in the ground state of harmonic 
chains similar to the one we studied here 1121 ■ 



VI. CONCLUSIONS 



FIG. 10: Dynamics of entanglement for oscillators in an initial 
separable state r = 2. The dynamics of non-resonant oscilla- 
tors along with the resonant case is shown for (a) ohmic and 
(b) sub-ohmic environments at T/fi — 10 with 70 — 0.1. 



to the action of the environment. However, entanglement 
decays much faster for non-resonant oscillators and the 
state becomes disentangled in a finite time (SD). We can 
also observe that the sub-ohmic environment can retain 
entanglement a bit longer than the ohmic one. This is 
due to the fact that the bare coupling (which is later 
cancelled by the coupling induced by the environment) 
produces more entanglement at short times in the sub- 
ohmic case. One can also notice little differences between 
the evolution corresponding to non-resonant oscillators 
with higher (or lower) frequencies. Indeed, this is due 

to the fact that the virtual interaction, c_| , depends on 

the square of the frequencies and not upon the differ- 
ence between them. As both virtual oscillators approach 
an equilibrium state (which is characterized also by non- 
vanishing correlations between them) the final entangle- 
ment turns out to be independent of the initial state. 
Thus, generic fate of entanglement at sufficiently high 
temperature is sudden death. 




T/n 



FIG. 11: Asymptotic entanglement between non-resonant os- 
cillators as a function of the temperature for different detun- 
ings A = — u!2 (ohmic spectral density). We can find a 
critical temperature for every detuning below which the state 
is entangled. They are compared with the entanglement for 
r = and A = 0. 



We presented a complete study of the evolution of the 
entanglement between two oscillators interacting with 
the same environment. We extended the analytical and 
numerical results previously presented in |24j . We con- 
sidered two related models for the interaction between 
the system and the environment: one where the coupling 
is through position and another where the coupling is 
symmetric in position and momentum. In both cases we 
used an exact master equation as our main analytical 
tool. For position coupling we presented a phase dia- 
gram valid for ohmic and sub-ohmic spectral densities, 
and we showed that it contains three phases (SD, NSD 
and SDR). For both spectral densities the phase diagram 
is qualitatively the same. The main difference is that the 
sub-ohmic environment tends to enhance the amplitude 
of the entanglement oscillations (which is due to the fact 
that the asymptotic state induced by a sub-ohmic envi- 
ronment has larger squeezing than the one corresponding 
to the ohmic case). On the other hand, we showed that a 
qualitatively different phase diagram emerges when the 
coupling is symmetric. In that case, the SDR phase is 
absent and the asymptotic entanglement does not oscil- 
late. Our results clearly show that initial separable states 
can get entangled and that initially entangled states can 
suffer from sudden death. 

For position coupling, we showed that there is a range 
of temperatures where SD never occurs. In fact, this is 
the case for T < Tq where Tq is the temperature where 
the position dispersion of the a;+ oscillator becomes iden- 
tical to the one corresponding to vacuum (below Tq such 
dispersion is smaller due to squeezing). On the other 
hand, for symmetric coupling the SD phase is present for 
every temperature. 

Our results can be extended in several ways. In fact, 
along the paper we have focused on the case where the 
renormalized oscillators do not interact, but our analysis 
can be applied to the cases where C12 does not vanish. 
In this ri+ 7^ ri_, the boundaries between dif- 

ferent phases change slightly but the phase diagram re- 
mains qualitatively the same. For example, the formula 
for Vcrit given in(22) tells us that when fl^ ^ f2_ r^rit 
may be non-zero even if the state of the oscillator is 
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not squeezed. In this case the resource for the asymptotic 
entanglement is, not surprisingly, supplied by the inter- 
action. On the other hand, when the coupling is sym- 
metric the asymptotic evolution does not change consid- 
erably by adding an interaction between the oscillators. 
Indeed, in this case, we always have Mfl = m^uj^ = mu. 
Therefore, r^rit vanishes and there are no entanglement 
oscillations in the long time regime. There is only one ex- 
ception to this rule: If one introduces a non-symmetric 
coupling between the oscillators, i.e. C12 7^ £12, then we 
get MQ. ^ m_Ci;_. Our results also change if the initial 
states of the system is mixed. However, the change in 
the phase diagram is simple to understand. In fact, the 
degree of purity of the initial state is characterized by the 
product 5x_5p^, which only enters in the expressions of 
Scrit and changes the mean value of the final entangle- 
ment (as seen in eqs. (19 1 and (23)). It is simple to see 



that the entanglement achieved for pure sates is grater 
than the one obtained for mixed initial states with the 
same degree of squeezing. The phase diagram for mixed 
states can be simply obtained from the one correspond- 
ing to pure states by shifting the curve Scrit to the right. 
This has the effect of moving upwards the horizontal axis 
(see Fig. [TJ . As a consequence, the value of Tq changes 
and the low temperature NSD island may disappear de- 
pending on the degree of impurity of the initial state. 

The existence of asymptotic entanglement between res- 
onant oscillators xi and X2 can be understood in terms 
of the following quantum-optical analogy: We can think 
of these oscillators as two modes of the electromagnetic 
field. The evolution of such modes, interacting with the 
environment is equivalent to the following sequence of 
operations: i) a 50/50 beam splitter supperpose both 
modes (creating the x± oscillators out of the original 
ones), ii) while one of the output modes (x_) evolves 
freely, the other is replaced by a new one with disper- 
sions along its quadratures given by the equilibrium val- 
ues (this operation entirely replaces the interaction be- 
tween and the environment), iii) another 50/50 beam 
splitter is applied (which gives rise to the final state of 
the Xi^2 oscillators out of the virtual x± ones. Follow- 
ing [33] we can conclude that the non-classicality at the 
output modes (after the second beam splitter) must arise 
from some form of non-classicality at the input. This can 
exist if the equilibrium state has some degree of squeez- 
ing (this is the case for position coupling) or if the initial 
state is non-classical (either entangled or squeezed). The 
(pure) initial state which is least favorable for produc- 



ing entanglement at the output are coherent states. The 
condition for the existence of entanglement in the final 
state for such initial states is r^^it > l/21n(2Aa::+Ap+). 
Thus, to fulfill this condition we need the environment 
to produce an equilibrium state where the variance of 
one of its quadratures is smaller than the vacuum limit, 
i.e. min{A2a;+, A^p+} < 1/2 (for m = l,ri_ = 1). We 
showed that this happens for position coupling and tem- 
peratures below Tq. Moreover, the above description of 
the problem enables us to draw a stronger conclusion: 
For symmetric coupling initial coherent states will never 
get entangled (even at intermediate times). This result, 
confirmed by our numerical simulations, can be seen as 
follows: For initial coherent states, the oscillators x+ and 
x_ are not squeezed initially. Moreover, x_ does not 
evolve for this type of coupling while the oscillator x+ 
will change its variances but will never become squeezed 
due to the nature of the interaction, which is symmetric 
in position and momentum. Therefore, the two modes 
will have vanishing squeezing during the entire evolution 
and, as a consequence the oscillators xi^2 will never be 
entangled. 

Finally, we studied the behavior of non-resonant oscil- 
lators by using both numerical and analytical tools. In 
this context we obtained a new master equation where 
the two virtual oscillators {x-^. and x_) are coupled. Both 
oscillators approach an equilibrium state where they are 
not only correlated but also may have slightly different 
variances. We showed the existence of an entangled state 
at very low temperatures. Thus, we conclude that the 
generic fate of entanglement in a finite temperature envi- 
ronment is not only to become independent of the initial 
state. Also, we showed that there is a low temperature 
threshold that depends on the detuning above which the 
entanglement undergoes a sudden death. This is prob- 
ably related with the entanglement studied in harmonic 
chains [H]. A more detailed analysis of the possible sce- 
narios for non-resonant oscillators will be presented else- 
where. 
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